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Deep XMM-Newton observation of a northern LMC field: 

I. Selected X-ray sources * 
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Abstract. First results from a deep XMM-Newton observation of a field in the Large Magellanic Cloud (LMC) near the north- 
ern rim of the supergiant shell LMC 4 are presented. Spectral and temporal analyses of a sample of selected X-ray sources 
yielded two new candidates for supernova remnants, a supersoft X-ray source and a likely high mass X-ray binary (HMXB) 
pulsar. From the fourteen brightest sources up to ten are active galactic nuclei in the background of the galaxy which can 
be used as probes for the interstellar medium in the LMC. From the three previously known HMXBs the Be/X-ray binary 
EXO 053 109-6609.2 was the brightest source in the field, allowing a more detailed analysis of its X-ray spectrum and pulse 
profile. During the pulse EXO 053109-6609.2 shows eclipses of the X-ray emitting areas with increased photo-electric absorp- 
tion before and after the eclipse. The detection of X-ray pulsations with a period of 69.2 s is confirmed for RX J0529. 8-6556 and 
a possible period of 272 s is discovered from XMMU J053011. 2-655122. The results are discussed with respect to individual 
sources as well as in the view of source population studies in the vicinity of the supergiant shell LMC 4. 
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1. Introduction 

Due to their proximity the Large and Small Magellanic Cloud 
(LMC, SMC) were always subject to X-ray surveys for satel- 
lites carrying imaging instruments. The most sensitive and 
most complete survey is available from ROSAT data obtained 
between 1990 and 1998 in the energy range 0.1-2.4 keV. In 
particular pointed observations with the PSPC detector with 
its large field of view covered in total ~59 and ~18 square 
degrees of th e LMC and SMC r egions on the sky (Haberl & 
Pietsch |1999| ; |Haberl et al.||2000| ). The ROSAT PSPC and HRI 
observations revealed about 1000 and 750 X-ray sources in the 
direc tion of the LMC and SMC, respectively (see also Sasaki 
et al. |2000a[ |bj). 

The large number of X-ray sources in the Magellanic Cloud 
fields allows population studies on rich samples of various 
kinds of source types powered by different mechanisms. An 
unusually high concentration of Be/X-ray binary systems was 
found in the SMC ( t-[aberl & Sasaki| |2000| ). Detectors sensi- 
tive to higher energies on board ASCA and XTE detected pul- 
sations from many new Be/X-ray transients during outburst 
dYokogawa et al.| |2002| ). The high sensitivity of the XMM- 
Newton instruments promises the detection of pulsations to 



much lower X-ray fluxes ( |Sasaki et al.||2001| ) usually observed 
from Be/X-ray binaries in their low state. 

To study the source population in the LMC to flux limits 
below 10 -14 erg cm~ 2 s _1 a deep —60 ks XMM-Newton ob- 
servation was performed as part of the telescope scientist guar- 
anteed time. The selected field, aimed at RA = 05 h 31 m 20 s and 
Dec = -65°57'38", was chosen because of its location on the 



rim of the supergiant shell (SGS) LMC 4 ( |Meaburr^l980[ ) and 
the simultaneous coverage of three known high mass X-ray bi- 
naries (HMXBs). These are : 1) The Be/X-ray binary pulsar 
EXO 053 109-6609.2 discovered in outburst duri ng EXOSAT 



observations of the LMCX-4 region in 1983 (Pakull et al 



1985). X-ray pulsations with a pulse period of 13.7 s were de- 



tected for the first time in ROSAT data by pennerl et al.| ( |1996| ). 
2) An outburst observed from RX J0529.8-6556by ROSAT and 
the detection of 69 s X-ray pulsations led to the suggestion 
by |Haberl et al.| ( |1997| ) that RX J05 29.8-6556 is also a Be/X- 
ray binary. Optical spectroscopy by Negueruela & Coe ( 2002} ) 
confirmed this. 3) The highly variable source RXJ0532. 5-6551 
proposed as first HMXB in the LMC powered by accretion 
from the wind of an OB supergiant companion ( |Haberl et al 



1995) which is supported by optical spectroscopy (Negueruela 



&Coe 2002) 
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* XMM-Newton is an ESA Science Mission with instruments and 
contributions directly funded by ESA Member states and the USA 
(NASA) 



Active galactic nuclei (AGN) behind nearby galaxies are of 
interest for two main reasons. They ca n be used to define a pre- 
cise reference coordinate system (e.g. Anguita et al. 200C ) and 
provide line of sigh t probes to study the interstellar medium 
in the galaxies (e.g. Kahabka et al. 2001 ; Haberl et al. 2001). 
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The new X-ray observatories Chandra and XMM-Newton will 
largely increase the number of background AGN as first re- 



sults from LMC observations demonstrate (Haberl et al 



2001 



Dobrzycki et al.||2002| ). 

In this article first results from the deep XMM-Newton 
observation are presented. The paper concentrates on spectral 
and temporal analyses of a selected sample of X-ray sources 
which includes the three known HMXBs and new candidates 
for supernova remnants (SNRs), a supersoft source (SSS), a 
likely Be/X-ray binary pulsar and ten candidates for back- 
ground AGN. 

2. The XMM-Newton observation 



XMM-Newton (Jansen et al. 2001) observed the field located 



in the northern part of the LMC during satellite revolution #152 
on Oct. 7, 2000 from 03:36 to 22:43 UT with a gap in the mid- 
dle of about 90 min due to ground station hand-over. A prelimi- 
nary analysis of the EPIC data, which did not apply proper cor- 
rections for telescope vign etting and point spread function, was 
reported in fr-Iaberl ( 2002 ). Here results of the analysis of a se- 
lected sample of detected X-ray sources are presented utilizing 
the data collected with the European Photon Imaging Cameras 
(EPICs). The cameras are b ased on MOS (EPI C-MOS1 and 
-MOS 2, |Turner et al] |200l| ) and pn (EPIC-pn, |Struder et af 
2001 ) CCD detectors and are mounted behind the three X-ray 
multi-mirror systems ( Aschenbach et al. 200C| ) and observe si- 
multaneously. All cameras were operated in Full-Frame read- 
out mode providing data over a field of view of ~13' radius 
and used the medium filter to block out optical light. The data 
were processed using the XMM-Newton analysis package SAS 
version 5.3.3 to produce the photon event files and binned data 
products like images, spectra and light curves. 

2. 1. The LMC field 

Images were created in the energy bands 0.3-1.0 keV (soft), 
1.0-2.0 keV (medium) and 2.0-7.5 keV (hard) and smoothed 
using intensity-dependent adaptive filtering. Detector-intrinsic 
particle induced background from an observation with closed 
filter wheel (on April 5, 2000) was subtracted before each im- 
age was corrected for telescope vignetting and exposure. A true 
colour image was produced by combining the three images and 
using red, green and blue colour scales to represent the intensi- 
ties in the soft, medium and hard energy bands, respectively. 
This RGB image is presented in Fig. [j]. According to their 
spectral energy distribution X-ray sources of different type are 
characterized by different colours. Sources predominantly ra- 
diating in the soft band appear in red while hard sources are 
blue. Nearby bright sources outside the field of view (an order 
of magnitude brighter than the brightest source inside the field) 
cause arc-like structures in the images due to photons which 
are single-reflected on the mirror shells. Arcs caused by the 
bright star AB Dor (at ^35' angular distance from the optical 
axis) can be seen in the NW part of the image. The supernova 
remnant N63A (an extended source ^27' off-axis) located East 
produces extended arcs on the left side of the image. Another 
arc from the pair of bright SNRs (N49/N49B, ~34') is visible 



on the edge of the image to the SW. To first order the energy 
dependence of the single-reflection efficiency is similar to that 
of "normal" double-reflected photons and the spectrum is lit- 
tle altered. Therefore, the arcs from the SNRs show the typical 
red-green colours of this class of sources with thermal spectra. 

3. Selected X-ray sources 

Source detection based on sliding window and maximum like- 
lihood methods available in the SAS package yielded about 
150 discrete X-ray sources in the field. A source list with X- 
ray properties and a statistical analysis will be presented in a 
follow-up paper. Here a sample of twenty sources is investi- 
gated which is selected on the basis of brightness, hardness 
and angular extent. To easily identify the sources in Fig. |l[ 
their locations are marked in a grey-scale broad-band image 
in Fig. ^. Some information useful for the identification of the 
sources is provided in Table [j]. Sixteen of them are located 
within 15" distance to ROSAT sources as presented in Haberl 
& Pietsch (|1999b and ISasaki et alJ (pOOOaj) and thirteen have 



a likely optical counterpart in the USNO A2.0 catalogue (one 
additional bright star is listed only in the HST guide star cata- 
logue). Angular distances between the X-ray and optical posi- 
tions are listed in the table (delta) together with the optical R 
and B magnitudes given in the USNO catalogue. 

Finding charts produced from the DSS2 (red) im age for 
some of the sources presented here are published in Haberj 
( 2002 ), however with preliminary X-ray source positions 
which can differ of the order of a few arc seconds. Among the 
optically brightest objects in the sample three are likely late- 
type foreground stars given their red colours. The identification 
of the X-ray sources with these stars is supported by the X-ray 
spectrum (see below). Three other objects with B brighter than 
15 mag are identified with known HMXBs and the fourth with 
XMMUJ05301 1.2-655122, a candidate HMXB proposed by 



Haberl (2002), based on the X-ray spectrum and the brightness 



of the optical counterpart which is similar to those of the known 
HMXBs in the field. |Haberl (2002) also proposed two new 
candidates for SNRs which can be easily identified in Fig. [j] 
from their angular extent and colours. For the brightest four- 
teen sources a sufficient number of counts (>500 in pn) were 
collected to allow spectral and temporal analyses. For sources 
with extreme X-ray spectrum a smaller number of counts al- 
lows a classification of the X-ray source. Six such sources were 
selected including a likely faint SSS, two SNRs, two more fore- 
ground stars and the known HMXB RX J0532.5-6551. 

3. 1. EPIC X-ray spectra 

Source spectra averaged over the whole observation were ex- 
tracted from MOS and pn data in circular regions around the 
source positions with typical radii of 16"-22". For the ex- 
tended candidate SNRs larger radii of 30" and 45" were used. 
Background was taken from nearby circular regions chosen to 
be at a comparable distance to the read-out node of the pn de- 
tector (RAWY coordinate) in order to avoid detector-intrinsic 
background variations. Redistribution matrices available with 
SAS 5.3 were used in combination with effective area files pro- 
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Fig. 1. EPIC-pn RGB color image produced from data in three energy bands (red: 0.3-1.0 keV, green: 1 .0-2.0 keV and blue: 2.0- 
7.5 keV). Each individual image was smoothed with an adaptive intensity filter. Instrumental background from an observation 
with closed filter wheel (from revolution 59) was subtracted before vignetting and exposure correction. The XMM-Newton 
pointing was aimed at RA = 05 h 31 m 20 s and Dec = -65°57'38" and covers about 13' in radius. North is to the top and West to 
the right. The brightest source, located in the South, is the Be/X-ray binary system EXO 053109-6609.2. 



duced by "arfgen" which include corrections for vignetting and 
point spread function losses. The set of three spectra (MOS1, 
MOS2 and pn) was analyzed using XSPEC by simultaneously 
fitting a common spectral model where the normalizations be- 
tween the instruments were independent parameters to allow 
for inter-calibration uncertainties. Errors and upper limits were 
derived for 90% confidence levels. Except for the brightest 
source (EXO 053 109-6609.2) which requires a more complex 
model, single component models attenuated by photo-electric 
absorption were fit to the spectra. EPIC spectra from eight 



representative sources from the sample are shown in Fig. [3[ 
Various source types can be clearly distinguished from their 
different energy distribution in the 0. 15-10.0 keV energy band. 

3.1 .1 . High Mass X-ray Binaries 

Although the Be/X-ray pulsar EXO 053 109-6609.2 is known 
since almost 20 years, little is known about its X-ray spectrum. 
The small angular separation to LMC X-4 made it impossible 
to take spectra with non-imaging instruments and only mea- 
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Table 1. Sample of selected sources detected in the EPIC observations. 
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Fig. 2. EPIC-pn 0.3-7.5 keV grey-scale image with the sources 
investigated in this work marked. 



surements with the EXOSAT Low Energy imaging detector 
with different filters at energies below 2 keV were available 
( Brunnei 1987). In the 2-20 keV band a spectrum of low sta- 



tistical quality was obtained from the coded mask SL2 XRT 
experiment which was consistent with a power-law with a pho- 



ton index 1.0^° 5 an d an absorbing hydrogen column density 
of 7±2. 6 xlO 22 cm' 2 ( Hanson et al.]|l989| ). 



The XMM-Newton observation yielded the first 0.2-12 
keV spectrum (Fig. |) of EXO 053 109-6609.2. The spectrum 
is complex and cannot be represented by the simple power- 
law (neither by including a high-energy cutoff) model sufficient 
to characterize most Be/X-ray binaries. Below 1 keV clearly 
a soft component is detected with features too narrow to be 
fit with continuum spectra as produced by e.g. blackbody or 
multi-temperature disk-blackbody emission. The only accept- 
able model was found to be thermal emission from hot plasma 
in collisional equilibrium (VMEKAL model in XSPEC, Mewe 
et al. 1985 ) with a temperature kT of 0. 1 keV and reduced abun- 
dances (0.29 solar) of elements heavier than oxygen. Above 
3 keV a highly absorbed (~10 23 cm~ 2 ) power-law gives an 
acceptable fit consistent with the SL2 results. However, both 
components described so far do not contribute significantly to 
the flux between 1 and 2 keV. Reducing the absorption in the 
power-law component yields too high flux around 2-3 ke V and 
does not fit. Several components were tested in the model to 
account for the emission between 1 and 2 keV. Neither black- 
body, disk-blackbody nor a second VMEKAL give acceptable 
results. 

The best fit was obtained by adding another power-law 
component which is attenuated at low energies by the same 
amount of absorption as the thermal plasma component. The 
power-law index was forced to be the same as the index of the 
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Fig. 3. EPIC spectra of selected X-ray sources in the LMC field. For each source the best fit to the three EPIC spectra is 
shown in the upper panels (pn is always the upper spectrum and the two MOS spectra overlap). The four HMXBs including 
the new candidate XMMU J05301 1.2-655122 are shown on the left. The right column illustrates the variety of X-ray spectra 
from sources of different nature. The spectra of candidate AGN show steeper power-laws than HMXBs. Spectra of the candidate 
SNR XMMU J053226.2-655352 and foreground stars exhibit little flux above 2 keV but can be distinguished by the amount of 
absorption. The softest spectrum is detected from the new candidate SSS XMMU J053056.2-654809 which can be represented 
by a blackbody model with kT~50 eV. For this case only the pn data yields sufficient statistics to accumulate a spectrum. 
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Table 2. Spectral fit results. 
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spectral fit to pn spectrum only 

the low absorption indicates a low energy excess which is probably caused by contamination from EXO 053 109-6609.2 
located in arc structure caused by single reflected photons from nearby bright SNR which probably causes contamination 
Spectral model; MC: multi-com ponen t (see text), PL: power-law, BB: blackbody, MK: thermal emission from plasma in collisional equilib- 



Mewe et al 



1985] ) 



rium (Mewe-Kaastra 

5 0.2-10 keV, determined from the pn spectrum except for EXO 053 109-6609.2 (MOS1) which is located on a CCD gap in the pn 

6 0.2-10 keV, intrinsic luminosity with Nh set to for sources located in the LMC assuming a distance of 50 kpc 



high absorption power-law. The best fit values for the photon 
index and the two column densities are listed in Table |. The 
normalization of the low-absorption component relative to the 
high-absorption component is 0.23. The two power-law com- 
ponents may originate from the same emission region partially 
covered by different amounts of absorbing matter. Partially 
may be interpreted here as spatially or temporally varying. 
The latter is supported by spectral variations seen during the 
pulse period (see below). Adding a narrow Fe emission line 
(6.32±0.06keV) with an equivalent width of 329 eV improved 
the x 2 by 13 to the best fit with x 2 =520 for 336 degrees of 
freedom. The model spectrum is illustrated in Fig. Q 

A soft spectral component is so far only detected from a few 
HMXBs, most of them located in the Magellanic Clouds (e.g. 
Kohno et alj |2000| ; |Yokogawa etat] |2000| ). In the Milky Way 
the HMXBs are concentrated in the galactic plane and the large 
interstellar absorption strongly attenuates the soft X-ray ener- 
gies. Therefore, the HMXBs in the Magellanic Clouds are ideal 
to investigate the soft part of the spectrum. No pulsations were 
detected in the soft spectral component seen from RX J0059.2- 
7138 with ASCA (Kohno et al. 2000) which suggests emission 
from a relatively large region. The emission measures derived 
for EXO 053 109-6609.2, if the thin plasma model is adopted, 




Fig. 4. Best fit model spectrum for EXO 053 109-6609.2. It 
comprises of two power-law components with equal slope ab- 
sorbed by different amounts of column density ("partial ab- 
sorber") and a thermal plasma emission component. For pa- 
rameters see Table 0. 



dKohno et al.||2000[ ) support this. A possible origin of the soft 
emission may be the disk-like dense stellar wind around the 



Be star which is illuminated by the strong X-ray source with 
luminosities of the order of 10 38 erg s _1 . However, to differen- 
tiate between emission lines produced in collisionally or photo- 
ionized plasma high resolution grating spectra are required. 
E.g. Chandra observations of the HMXB CenX-3, a close sys- 
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tem with supergiant primary, revealed emission lines excited 
by both kinds of processes ( Wojdowski et al. 2002). 

From RX J0532.5-655 1 and RX J0529.8-6556 only ROS AT 
PSPC spectra in the energy band 0.1-2.4 keV were avail- 
able. The spectra are of low statistical quality and the nar- 
row energy coverage did not allow to determine power- 
law in dex a n d col u mn den si ty wi th high accuracy (Haberl 
et al. |!995j |1997| ). |Uaberl| ( |2002| ) proposed a new candi- 
date HMXB based on its hard X-ray spectrum and the pres- 
ence of a likely optical counterpart on the DSS2 (R) im- 
age with similar brightness as the counterparts of the other 
two nearby Be/X-ray binaries. Detection of a possible pulse 
period (see below) for XMMUJ05301 1.2-655122 leaves lit- 
tle doubt about the HMXB and most likely Be/X-ray binary 
nature of this source. Fitting an absorbed power-law to the 
spectra of XMMUJ05301 1.2-655122, RXJ0529. 8-6556 and 
RXJ0532.5-6551 yields acceptable fits (reduced x 2 0.9-1.3) 
with power-law indices in the range 1.0-1.4 (Table The 
EPIC spectra of these sources are shown together with the best 
fit model in Fig. |. X MMU J05301 1.2-655122 was also de 



tected in ROSAT data (Haberl & Pietsch 1999; Sasaki et al 



2000a| ) with a count rate corresponding to a flux which is a 
factor of 2.2 below the XMM-Newton measurement (assuming 
spectral parameters derived from the EPIC spectra). 



3.1.2. AGN 

The majority of the X-ray brightest sources in the field (ten out 
of fourteen) most likely consists of AGN in the background of 
the LMC. Their X-ray spectra can be represented by power- 
laws with indices between 1.8 and 2.3 with only one case with 
index 1.4 (Table |J), typically for this source class. Moreover 
they do not possess an optical counterpart with B brighter than 
18th magnitude in the USNO A2.0 catalogue (Table [[]). As 
an example the EPIC spectra of XMMUJ053150.9-655616to- 
gether with the best fit model are plotted in Fig. |5| 



3.1 .3. Supernova remnants 



Two new SNRs were proposed by frlaberl (2002) from their 
spatial extent and X-ray colours (Fig. |lj) in the EPIC data. 
Both SNRs were detected in the ROSAT PSPC data (sources 
190 and 205 in frlaberl & Pietsch] |1999| ) with indication for 
spatial extent. However, due to the low statistical quality of 
the data the sources were not classified as SNR candidates. 
XMMUJ053226.2-655352 is sufficiently bright to allow the 
investigation of the X-ray spectrum which shows little flux 
above 2 keV. Fitting a MEKAL model results in a best fit 
temperature kT = 0.68 keV consistent with radiation from 
a thermal SNR. For the Fe abundance a best fit value of 
0.4 l+o ii was derived. The best fit model to the EPIC spectra of 
XMMU J053226.2-655352 is shown in Fig. | Due to the lower 
number of detected photons from XMMUJ053034.9-655653 
the temperature was fixed at 0.7 keV in the fit to the EPIC spec- 
tra. Also in this case a low Fe abundance was deduced (upper 
limit of 0.06). 



3.1.4. Foreground stars 

The EPIC spectra of sources identified with foreground stars 
also do not show significant emission above 2 keV. However 
the absorbing column density is lower than that of sources lo- 
cated in the LMC clearly distinguishing them from SNRs (in 
addition they are unresolved X-ray sources). Spectra of three 
foreground stars were investigated and the results are listed in 
Table ||. All three spectra can be represented by emission from 
hot thermal emission with temperature between 0.5 and 0.7 
keV and low Fe abundance of 0.1-0.2 solar. For the absorption 
column densities only upper limits could be derived consistent 
with nearby objects. The X-ray spectra suggest emission from 
late-type stars with active corona which is also consistent with 
their B-R index derived from the magnitudes listed in USNO 
A2.0. As an example the spectra of XMMU J052924.3-655724 
are plotted in Fig. 0. 



3.1 .5. Supersoft source 

The faint source XMMU J053056.2-654809 is characterized by 
red colour in the image of Fig. [j] indicating an extremely soft 
X-ray spectrum. Only the pn detector is sensitive enough to 
allow accumulation of an energy spectrum which is only de- 
tected below 0.5 keV (see Fig. ||). A formal fit with a black- 
body model yields a temperature of 49 eV which is typical for 
supersoft sources, AM Her systems, isolated white dwarfs or 
isolated X-ray dim neutron stars. Its luminosity of 1.3 xlO 34 
erg s _1 for a distance of 50 kpc is too high to originate from 
an AM Her type cataclysmic variable in the LMC (for lumi- 
nosities ofJhe_s^ft_blackbody component in AM Her systems 
see e.g. |Ramsay eTaI||l994| ). XMMU J053056.2-654809 is not 
contained in the ROSAT catalogues of the LMC and there- 
fore archival PSPC data (45 ks observation with ID 200692P) 
was investigated and an upper limit of 7.9 x 10~ 4 cts s _1 (3a, 
0.1-0.6 keV) was determined. Folding the spectral model fit- 
ting the EPIC pn data through the PSPC response an expected 
count rate of 2.6 xlO -2 cts s _1 is obtained which implies that 
XMMUJ053056.2-654809 was at least a factor of 33 fainter 
during the ROSAT observation in October/November 1991. 
The high variability excludes an explanation as dim isolated 
neutron star or single white dwarf in the foreground. An AM 
Her system located in the Milky Way is expected to be brighter 
in the optical than the close USNO catalogue object. Therefore, 
XMMUJ053056.2-654809 is considered most likely to be a 
variable SSS in the LMC. 



3.2. X-ray pulsars 

Temporal analysis was performed on the data of the HMXBs 
and the other eleven brightest sources in the sample. Fast 
Fourier power spectra were produced from light curves in the 
0.3-7.5 keV energy band. A folding analysis was then per- 
formed around frequency peaks identified in the power spectra. 
Pulsations were detected only from the known HMXB pulsars 
and from the new candidate HMXB which is discussed in the 
following. 
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3.2.1. EXO 0531 09-6609.2 

For EXO 053 109-6609.2 a mean pulse period of 13.66817(1) 
s was derived (Ict error given for the last digit). The power 
spectrum exhibits also peaks at the first and third harmonic, 
consistent with the highly non-sinusoidal pulse profile (Fig. ||). 
Pulsations are seen above and below 4 keV confirming that 
the two spectral power-law components (Fig. [|) originate both 
from EXO 053 109-6609.2. Due to the strong re-distribution of 
events to lower energy-channels at energies below ~1 keV in 
the EPIC pn CCD it is not possible to verify if the soft thermal 
plasma emission below 0.5 keV shows pulsations or not. For 
a spectrum steeply falling towards lower energies the contribu- 
tion from photons above 0.5 keV dominates the events below 
channel 100 (nominally corresponding to energy 0.5 keV) and 
therefore pulsations seen in the flux above 0.5 keV will also be 
seen in the lowest channels. 

To investigate possible spectral changes during the pulse 
period, hardness ratio plots were produced by creating light 
curves folded on the pulse period in different energy bands 
and dividing the obtained count rates as function of pulse 
phase. The folded light curves from energy bands 0. 15-0.4 keV 
and 0.4-2.0 keV are plotted in Fig. ^ (left) together with the 
hardness ratio. The hardness ratio shows remarkable increases 
around the two intensity dips which resemble the behaviour of 
an X-ray source viewed through a dense atmosphere before and 
after it is eclipsed. The facts that the eclipse seen here is not to- 
tal and that two eclipses are seen about 0.5 in phase apart, sug- 
gest that we see two emission regions in EXO 053 109-6609.2, 
most likely from the accreting magnetic poles of the neutron 
star. Before the pole is self-eclipsed by the neutron star or the 
accretion column the emission is viewed through dense matter, 
most likely matter from the accretion column above the pole 
which attenuates the emission more strongly at lower energies. 
During eclipse (when the emission is probably dominated by 
the other pole) the hardness ratio remains constant within the 
errors. Hardness ratios which involve emission from energies 
above 2 keV (which is less affected by absorption and more 
sensitive to changes in the intrinsic spectrum) are plotted in 
Fig. H (right). They show that the spectrum softens during the 
eclipses. This can be explained by a spatial temperature vari- 
ation in the emission region with the temperature increasing 
with height above the neutron star surface. The cooler regions 
become only visible during closest approach of the magnetic 
axis to the line of sight. Such a picture is expected from mat- 
ter in the accretion column shocked at a certain height above 
the neutron star surface and cooling during the final settlement 
onto the neutron star. 

Comparison of the pulse profiles of the XMM-Newton ob- 
servation with those obtained from BeppoSAX data (Burderi 
et al. 1998) shows a major difference which is best visible in 
the 2-10 keV band. During the BeppoSAX observation the pro- 
file shows a strong peak after the deeper eclipse while after the 
secondary eclipse the intensity remains as low as during the 
eclipse. During the XMM-Newton observation this behaviour 
is mirrored around the secondary eclipse (at phase ^0.55 in 
Fig. ||) with low intensity before and high intensity after the 
eclipse. This behaviour is probably related to the overall source 



intensity. The pulse profile obtained during a long ROSAT ob- 
servation ( Dennerl et al. 1996]) when the source was in a low- 



intensity state (1 x 10 erg s ) resembles those of the XMM- 
Newton observation (with observed luminosity of 1.7 xlO 35 
erg s _1 in the 0. 1-2.4 keV band). During the BeppoSAX obser- 
vation the source was in outburst with an observed luminosity 
in excess of 10 37 erg s _1 compared to 2.1 x 10 36 erg s _1 mea- 
sured by EPIC (note the luminosity given in Table || is corrected 
for absorption). It should be noted that no attempt was made in 
the presented analysis to determine secular period changes and 
the pulse profiles based on a mean pulse period may be smeared 
out somewhat. 



3.2.2. RXJ0529.8-6556 and 

XMMUJ05301 1.2-6551 22 

The pulse period of 69.232(2) s detected in the EPIC data 
from RXJ0529.8-6556 confirms the discovery by ROSAT. 
For the second brightest X-ray source in the field during 
the XMM-Newton observations, XMMUJ05301 1.2-655122, 
probable pulsations were found with a period of 271.97(5) 
s. From the folding analysis the \ 2 test yields a formal sig- 
nificance of 4.6 a for the probability that this detection is 
real. Together with the X-ray spectrum and the brightness of 
the proposed optical counterpart the pulsations strongly sup- 
port the identification of XMMUJ05301 1.2-655122 as fourth 
HMXB system in the observed field. Power spectra and light 
curves folded with the pulse period for RXJ0529. 8-6556 and 
XMMU J05301 1.2-655122 are shown in Fig. | 

4. Source populations 

Chu ( |1997[ ) used a sample of SNRs in OB associations in the 
LMC to study how interstellar environments affect the physi- 
cal properties of SNRs. She found that SNRs in HII regions 
all show the three classical SNR signatures (bright X-ray emis- 
sion, non-thermal radio emission and enhanced [S II]/Ha ratio) 
while some SNRs in super-bubbles show only X-ray emission. 
The observed LMC field presented here is located at the north- 
ern rim of the SGS LMC 4 and one reason that the two new 
X-ray selected SNRs were not detected yet at radio or opti- 
cal wavelengths may be the explosion of their progenitors into 
the hot and low density interior of the SGS. Very likely the 
formation of the progenitors was triggered by the expanding 
shell and the SNRs are now shocking the inner walls of the 
shell (|Chu||l997[). A relatively large number of SNRs with X- 



ray luminosities around 10 34 erg cm~ 2 s _1 may therefore still 
be detectable in XMM-Newton observations of the Magellanic 
Clouds and in particular inside the supergiant shells in the 
LMC. Promising candidates may be found among the sources 
listed in the ROSAT catalogues (|Haberl & Pietsch||1999| ; Haberl 
et al. 2000; Sasaki et al. 2000a with indication for spatial 
extent but without sufficient number of detected photons (and 
therefore a low likelihood for the extent) to allow a classifica- 
tion of the sources. 

The SMC is peculiar in its high number of Be/X-ray bi- 
naries (both in absolute terms given the size of the galaxy and 
relative to the number of supergiant HMXB s and SNRs, Haberl 
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Fig. 5. Power spectra and pulse profiles of HMXB pulsars produced from EPIC-pn data in the 0.3-7.5 keV energy band. 



& Sasaki 2000; Yokogawa et al 2002) which may indicate an 
increased star forming rate about 10 Y years ago (the time which 
elapses between the f ormation of a massive b inary and its evo- 
lution into a HMXB, van den Heuvej 1983 ). The star forma- 
tion must have largely declined since then as the lower number 
of descendants of young mas sive stars like supergian t HMXBs 
and type-II SNRs indicates ( Yokogawa et al.| 2002 ). With re- 
spect to the space density of HMXBs the observed XMM- 
Newton field and its close neighborhood in the northern area of 
the SGS LMC 4 are similar to that of the SMC. However, the 
existence of a relatively high number of younger SNRs and the 
supergiant HMXB system RXJ0532. 5-6551 is consistent with 



a more constant star formation rate over the last 10 7 years in 
this part of the LMC 4 region. The stochastic self-propagating 
star formation (SSPSF) scen ario proposed to expla in the pro- 
duction of supergiant shells ( ^eitzinger et al. 1981 ) predicts a 
gradient in the age of stars, from young stars at the edges to 
older ones inside. However, no such gradient was found for 
stars in the interior of LM C 4 with stellar ag es in a narrow 
range of 9 Myr to 16 Myr ( Braun et al. 1997 ) leading to the 
conclusion that LMC 4 may have formed without contribution 
from SSPSF, although the ring of young associations and H II 
regions around the edge were triggered by the events inside 
the SGS. The X-ray source census of HMXBs and SNRs in the 
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Fig. 6. Folded light curves obtained from EPIC pn data from EXO 053109-6609.2 in different energy bands together with the 
hardness ratios (count rates in hard band / count rates in soft band). 



LMC 4 region is fully consistent with this picture and the Be/X- 
ray binaries may well belong to the stellar population found by 



Braun et al. (1997). 



The majority of X-ray sources in the observed field are 
most likely AGN behind the LMC. Their spectra are charac- 
terized by power-laws attenuated by photo-electric absorption 
due to interstellar gas in the Milky Way and the LMC along 
the line of sight. From HI 21 -cm maps the Galactic contribu- 



tion is measured to 5.7 x 10 cm ( Dickey & Lockman 199C ) 



while the LMC col umn density var ies between 5x10 cm 
and 9xl0 20 cm~ 2 ( Kim et al. 1998 ) across the field. From the 
ten AGN candidates about half have X-ray measured column 
densities well above the HI values. Since the absorption in 
the X-ray band is dominated by the metals in the absorbing 
gas lower metallicities adopted in the spectral fits (0.5 solar for 
the LMC, Russell & Dopita 1992) would lead to correspond- 
ingly higher hydrogen column densities. Similar results, but 
with larger excess column densities, were obtained for AGN 
detected in the 30 Doradus region and Haberl et al. (2001 ) con- 
cluded for the 30 Doradus region that either large amounts of 
H II and/or H2 or higher metal abundances exist along the line 
of sights through the LMC. The case of the XMM-Newton field 
presented he re is less clear. No dense clouds are visible on Ha 
image s (e.g . Domgorgen et al. 1995 ) or CO maps (Yamaguchi 
et al. 2001 ) in the observed field at the northern rim of the 
LMC 4 SGS and in the sample of candidate AGN the excess 
absorption may be explained by AGN intrinsic obscuration. 
Indeed relatively large fractions of obscured AGN were found 



in deep surveys (e.g. fVIainieri et al. 2002), however, mainly 
among the X-ray and optically fainter ones. The ten brightest 
AGN in the LMC field have 2.0-10.0 keV fluxes higher than 
3xl0 -14 erg cm -2 s _1 and six of them likely counterparts 
with R magnitudes (from USNO A2.0) between 17 and 18, 
confining them to th e upper left par t of th e flux - magnitude 
diagram of Fig. 6 in Mainieri et al. ( 2002 ). There objects are 
located which show a much lower fraction of obscured AGN. 
To help to identify those unobscured AGN which can be used 
to probe the interstellar medium in the LMC optical spectra are 
required. 

5. Summary 

The deep XMM-Newton observation of a northern field in the 
LMC allows to uniquely classify the detected X-ray sources to 
flux levels of ^10~ 14 erg cm~ 2 s _1 from their X-ray proper- 
ties alone. X-ray spectra obtained by the EPIC cameras over 
a relatively broad energy band (0. 15-12.0 keV) emphasize the 
various energy distributions observed from the different types 
of X-ray sources. A first analysis of twenty selected sources, 
presented here, yielded the discovery of a previously unde- 
tected supersoft source, revealed two ROSAT sources with in- 
dication for spatial extent as supernova remnants and rendered 
the detection of a possible pulsation in the X-ray flux of a 
candidate HMXB in the LMC. This increases the number of 
known HMXBs in the observed field to four. Among the four- 
teen brightest sources three HMXBs, one foreground star and 
up to ten AGN were found. Particular results on individual ob- 
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jects are discussed in Sect. 3 while in Sect 4 the analyzed sam- 
ple of X-ray sources is considered in the view of source popu- 
lation studies in the Magellanic Clouds and nearby galaxies in 
general. 
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